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Coherent spin rotation in the presence of a phonon-bottleneck 
effect 
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Abstract. - A characteristic of spin reversal in the presence of phonon-bottleneck is the deviation 
of the magnetization cycle from a reversible function into an opened hysterezis cycle. In recent 
experiments on molecular magnets (e.g. V15 and RU2), the zero-field level repulsion was sufficiently 
large to ensure an otherwise adiabatic passage through zero-field and the magnetization curves can 
be described by using only a phonon-bottleneck model. Here, we generalize the phonon-bottleneck 
model into a model able to blend the non-adiabatic dynamics of spins with the presence of a non- 
equilibrium phonon bath. In this simple phenomenological model, Bloch equations are written 
in the eigenbasis of the effective spin Hamiltonian, considered to be a two-level system at low 
temperatures. The relaxation term is given by the phonon-bottleneck mechanism. To the expense 
of calculus time, the method can be generalized to multi-level systems, where the notion of Bloch 
sphere does not apply but the density matrix formalism is still applicable. 
EPL, 87 (2009) 57010; 23 September 2009; doi: 10.1209/0295-5075/87/57010 



Introduction. — Molecular magnets [1] or less com- 
plex systems containing diluted spins [2] have been ex- 
plored intensively in the recent years for their potential 
application in information technology. In addition, despite 
their macroscopic dimensions, crystals containing mag- 
netic molecules show measurable manifestations of fun- 
damental quantum mechanical phenomena at large scale. 
This is due only to the fact that the magnetic molecules are 
relatively well separated from each other and their quan- 
tum properties are amplified by the large number con- 
tained in a tiny monocrystal. Consequently, they show 
remarkable phenomena, like Quantum Tunneling of the 
Magnetization (QTM) measured first [3] in spin S = 10 
systems (Mni2, Fes) and Berry Phase quantum interfer- 
ence [4] (Fes). This large spin molecular magnets show sig- 
nificant anisotropy barriers protruded by tunneling chan- 
nels which are generated by a transverse magnetic field or 
anisotropy terms. A sweeping longitudinal field (that is, 
parallel with the main spin quantification axis) will force 
the spin to sequentially visit the tunneling channels and 
obey the Landau-Zener tunneling mechanism each time. 
Consequently, distinct jumps are generated in the hys- 



terezis cycle [3] allowing a direct measure of the tunneling 
gap [5]. 

Another type of molecular magnets consist of species 
having a small total spin, e.g. 5 = 1/2 (V15 [6], Vq [7]) 
or S = 1 (Feio) [8]. For a two-level system (S = 1/2) 
one can have S+, S- transverse terms generated by the 
complex structure of the molecules (like Dzyaloshinski- 
Moriya interactions in V15 [9]) leading to a level repulsion 
in zero field A (similar to a tunneling channel in large spin 
molecules). Due to the low value of the spin, there is no or 
low anisotropy barrier between states with opposite pro- 
jections of Sz- For A large enough to allow the spin to flip 
adiabatically under swept field, the magnetization cycle is 
reversible at thermal equilibrium. However, experiments 
performed in V15 and other molecules [9, 10] have shown 
an opening of the hysterezis cycle despite the absence of 
a barrier. Such phenomena have been interpreted within 
the frame of the phonon-bottleneck model which predicts 
a hysterezis of spin temperature in the presence of a fast 
sweeping field. For very small values of the level repulsion 
A, the phonon-bottleneck and the Landau-Zener tunnel- 
ing phenomena have to be treated together, which is the 
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main aspect of this article. 

Molecular systems are tunable, identical, two- or multi- 
level systems that can be produced in large numbers and 
therefore are potential candidates for qubit implementa- 
tion in quantum computing algorithms [11]. In practice 
however, spin qubits are very sensitive to decoherence 
mechanisms, mainly due to hyperfine and dipolar inter- 
actions with the background spins, and also to environ- 
mental thermal effects (phonon bath). Advanced chem- 
istry techniques allow synthesizing single crystals with a 
reduced amount of background magnetic elements leading 
to improved coherence times [2,12], whereas low tempera- 
tures reduce the phonon bath effect in usual experiments. 
In this paper we analyze the possibility of quenching ther- 
mal decoherence effects by operating the spins in phonon 
bottleneck (PB) conditions. 

Modeling. — The phonon-bottlcneck phenomenon is 
generated by an imperfect thermalization of a sample with 
the surrounding experimental setup. At low tempera- 
tures, the number of available phonons is strongly sup- 
pressed and the inherent limited sample thermalization 
with the cryogenic bath will delay the return to equilib- 
rium of sample's phonon bath. In the case of a magnetic 
sample, phonons are absorbed and re-emitted resonantly 
by the spin system. The timing of such a process, in the 
presence of a swept magnetic field, generates a bottleneck 
(strong absorption) followed by a phonon avalanche (de- 
layed strong re-emission) [13]. 

At low temperatures, phonons have a very low heat ca- 
pacity [13] ensuring that the lattice and the spins are at 
the same temperature Tg, which can be different however 
from the bath (cryostat) temperature T. As shown below, 
cycling of an external magnetic field can induce phonon 
bottleneck and a delay (hysteresis) in the dynamics of the 
spin temperature, which appears experimentally as a hys- 
terezis in the magnetic moment of a single crystal [9, 10]. 
In Ref. [13] a PB model based on the detailed phonon bal- 
ance is given and its implementation for the case of molec- 
ular magnets is presented in [9, 10]. In the case of a two- 
level spin system (e.g. S = 1/2), the model requires two 
fit parameters, a and A . The first one is a phenomeno- 
logical parameter related, among other, to sample's ther- 
malization and the second one is the level repulsion (or 
the tunneling gap) in zero field. Other theoretical models 
have been discussed as well in the literature [14]; they are 
able to simulate magnetic hysterezis but are less adapted 
to study the relation between sample thermalization and a 
phonon bath out of equilibrium. Also, recent microscopic 
modeling solving numerically the density matrix equation 
does give good quantitative description of experimental 
results [15]. 

A more intuitive model is discussed below and it con- 
sists in a modification of the well-known Bloch equations 
used in spin dynamics, to account phenomenologically for 
both the Landau-Zener process (not accounted for in the 
original model [13]) and the PB effect. 



Usually, a two-level system is described by the effec- 
tive Hamiltonian H e ff = —sa z /2 — Ao~ x /2 where a z , x 
are the Pauli spin matrices, A is the level repulsion and 
e = gpsH is the spin bias (ps is the Bohr magneton and 
H is the external magnetic field, always || to the z axis of 
H e ff). The resulting energy separation can be written as 
Ah = \J (A 2 + e 2 ) which is also the energy hv = Ah of 
phonons responsible for spin thermalization. In the pres- 
ence of a phonon bath in equilibrium, the spin dynam- 
ics can be described by the well-known phenomenologi- 
cal Bloch equations, particularly when a spin resonance is 
driven by a microwave excitation. The Bloch equations 
are equally useful when studying the effect of a ramping 
field on the dynamics of an isolated spin. For instance, 
when sweeping non-adiabatically e from — oo to +oo, the 
final value of (S z ) obeys the Landau-Zener [16] theory: 

Plz = 1 - cx P(- 2ftdf/dt )' wnere p lz is the probability 
for (S z ) to remain in the original state (up or down). 

In analogy to the classical case of magnetic resonance 
(A = 0), we will discuss first the Bloch equations in the 
diagonal basis of H e f / and how they can include the effect 
of phonon bottleneck. 

The magnetization components are given by Mi — 
(<7i) = Tr(pui) 1 i — x,y,z where o~i are the Pauli ma- 
trices and pis the density matrix. For a two level system 
one gets: 
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The magnetization dynamics is governed by Liouville 
equation dp/dt = —i/K[H(t), p] which can be written as 
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By evaluating pu ±p 2 2 and pi2±p2i one obtains the Bloch 
equations (with the addition of terms depending on A in 
absence of any relaxation or dephasing terms): 
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(2) 



Energy and phase relaxation are described phenomenolog- 
ically by additional decay terms, characterized by times 
T\ and T2 respectively. To this purpose, one has first to 
rotate Equations (2) into the eigenbasis of the spin Hamil- 
tonian to properly add the energy relaxation term. The 
level repulsion term Ao~ x /2 is equivalent with an applied 
field along x-axis which combines with the applied bias e 
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along z-axis to give the total magnetic field and the quan- 
tification axis of the cigcnbasis. It is therefore necessary 
to rotate the xyz coordinates around the y-axis by an an- 
gle = arctan(A/e). In the new x'y'z' coordinates, the 
averaged spin components (or magnetization) are: 



M x , = M x cos 9- M z sin 9 

My, = My 

M z , = M z cos 9 - M x sin 9 



(3) 



Time derivatives of the above set and Equation (2) lead 
to a new set of motion equations 
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One can introduce in Equations (4) the phcnomenological 
energy relaxation term along z'-axis and a dephasing term 
in the x'y' plane. The two processes are assumed to induce 
exponential decays, characterized by the relaxation time 
Ti and the dephasing time T 2 , respectively: 
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(5) 



where M e J is the magnetization value at thermal equi- 
librium. This set of equations, and in particular the en- 
ergy relaxation term, would correctly describe a typical 
spin-phonon interaction with the phonon bath at equi- 
librium. To describe correctly the particular case of a 
phonon-bottleneck driven process, the Ti term needs to 
be modified, as discussed below. 

Even at very low temperatures, the magnetic molecules 
are still coupled with the thermal heat bath and the spin 
system will relax towards equilibrium by exchanging en- 
ergy with the phonons. Only the phonons having the same 
frequency as the resonant frequency of the spins hu = Ah 
can do the energy transfer. The situation is schematically 
represented in Fig. 1. The spin system has a character- 
istic temperature T s defined as ni/n 2 = cxp(A# /fceTg) 
with ni,2 the out of the equilibrium occupation numbers 
of the ground and excited spin level respectively. The 
phonon bath has a temperature T p h defined by the Bose- 
Einstcin distribution pt pH £ = l/( ex P( £ /^-B^p?i) — 1). The 
bath temperature T represents the phonon temperature 
at equilibrium. The number of resonant phonons in an 
energy window Aid is [13] 
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where a(uj)dio — 3aj 2 duj/ (2ir 2 i> 3 ) is the number of phonon 
modes between u> and u> + du> per volume unit, v is the 
phonon velocity and Aui is the transition linewidth. As an 
example, let us consider the case of V15 molecule [6]. Tak- 
ing the typical values v=3000 m/s, T—0.1 K and Aw «500 
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Fig. 1: Illustration of the spin-phonon-bath coupling mech- 
anism: spins at temperature T s and heat capacity Ch are 
coupled to phonons at temperature T p h and heat capacity 
Cph « Ch, the energy transfer taking place with a rate 1/ti. 
Phonons are coupled to the cryostat at fixed temperature T, 
their relaxation rate being l/r p h- At very low temperatures, 
T s = Tph and the spin-phonon system relaxes toward T with a 
rate Cp h /{C H Tph) (see text). 

MHz, we find n T of the order of about 1CT 6 to 1(T 8 
phonons per molecules which is quite low. These phonons 
are rapidly absorbed by spins such as T s = T p h and the 
coupled spin-phonon system relaxes towards the bath tem- 
perature T. At low temperatures, absorption of resonant 
phonons generates a hole in phonon distribution between u 
and cu + dui, phenomenon known as phonon-bottleneck [17] 
and seen first experimentally in paramagnetic salts [18]. In 
such cases both the anharmonic phonon-phonon coupling 
and the phonon thcrmalization with the cryostat arc in- 
sufficient to prevent the hole formation in the phonon dis- 
tribution. Spin relaxation implies the transfer of energy 
from the spin system to the bath via phonons and can 
be treated simply by means of detailed balance equations. 
The spin population and the phonon number are given 
quantitatively by the following two differential equations 
describing the energy relaxation process: 
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where r p h is the relaxation time of phonons, Pi2,2i are 
the transition probabilities between the two spin levels 
and are directly proportional to the available number of 
phonons(cx px h and p^ h + 1 respectively). 
By using the notation 



1/Tl = (Pl2 + P21. 



eg 



(8) 



for the spin relaxation rate at equilibrium and the vari- 
able change (equilibrium quantities are indicated by the 
subscript cq) 



m - n 2 



-, y 



PTph - PT 



nieq - n 2eq pr - 1/2 

the balance equations can be rewritten as [13] 

dx/dt = (1 — x — xy)/ri 
dy/dt = —y/Tph + bdx/dt 

Here 6=5^ 



(9) 



(10) 



= "i+"2 tanh 2( Ray j th C and (J 

spin and phonon heat capacity respectively [13]. Solutions 
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of the system (10) are discussed in [19] for different values 
of T p h and b. Generally, in the vicinity of equilibrium 
(x — y 1, y — > 0), the spin relaxation is described by an 
exponential with a characteristic time 77, = T\ + (b+l)T ph . 
Since at very low temperatures Cjj ^ C p h and 6 > 1, 77, 
reduces to t h = br ph . For the very same reason, T s w T p/l 
[13], such that y = - — 1 Therefore, the second equation 
in (10) becomes 

<te = ll_-x 
at br p h x 

representing the relaxation law of the spin system in the 
presence of phonon-bottleneck phenomena, in the vicinity 
of equilibrium. The relaxation time tr depends on the 
applied magnetic field via the A#(e) dependence and is 
given by 

a . tan h 2 (-^-) (12) 



TH 



A 2 

H 



with a = 2tt 2 % ^ 3 A^r p / l /(3Aw) where N is the spin den- 
sity. 

Equation (11) allows modifying accordingly the phe- 
nomenological relaxation term in the Bloch equations (5) 
which was valid only for an equilibrated phonon bath. In 
the case of phonon bottleneck, eqs. (5) can than be rewrit- 
ten as 



= HT M v' ~ M * 
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Results and Discussion. — The modified Bloch 
equations (13) and eq. (12) allow simulating the spin mo- 
tion during a field sweep in the presence both phonon- 
bottleneck and Landau-Zener mechanisms (PB-LZ). The 
obtained hysterezis cycles do compare very well with the 
experimental ones in the case of large A (that is, without 
Landau-Zener non-adiabatic transitions), similarly to the 
PB-only method used in Refs. [6,9, 10]. 

To evaluate the effect of the Landau-Zener mechanism, 
magnetization cycles are calculated for values of A cor- 
responding to non-adiabatic and up to almost adiabatic 
conditions. Figure 2a shows (S z ) during half-cycles (field 
up-ramping) with the spin initially aligned to a negative 
field. The temperature is 0.3 K, magnetic field ramp rate 
is 0.1 T/s, attainable in current experimental conditions, 
a=l sK 2 and A=10~ 6 K and 10~ 5 K. Such values for 
a simulate an intermediate sample thermalization (com- 
pared with a ~l-3 sK 2 in the experiments performed on 
RU2 [10] or 0.09/130 sK 2 for good/bad thermalization re- 
spectively, in the case of the V15 experiments [9]). The 
dephasing time T2 is large (0.1 ms), so the simulations can 
evidentiate the spin coherent dynamics for a field interval 
up to 10" 5 T. 

The curves in fig. 2a give (S z ) for negative fields with 
values in between the zero temperature case and the Boltz- 
mann equilibrium case as a consequence of the phonon- 
bottleneck effect (0 < T s = T ph < T). The effect is even 
more visible after passing through zero field, and, for this 
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Fig. 2: (a) Half hysterezis cycles (field ramped up) for two 
values of A corresponding to almost adiabatic (thick line) and 
non-adiabatic (thin line) spin reversal. Due to the PB effect the 
spin is far from thermal equilibrium, (b) and (c) Spin dynamics 
in the two situations: adiabatic reversal almost contained in the 
xz plane and large spin oscillations due to the LZ effect. 



particular parameters, a valley of negative differential sus- 
ceptibility is seen (as in the experiments [10]). This is due 
to an increase in phonon energies as a result of spin de- 
excitations generating phonons with larger Tiuj than those 
absorbed in lower fields {phonon avalanche) . The last part 
of the curves shows a return to equilibrium. The PB-LZ 
method (14) also gives access to the spin trajectory, as 
shown in Fig. 2(b) and (c). In the adiabatic case (A=10 -5 
K) the spin trajectory is indeed approximately contained 
in the xz plane, defined by the field A and e. For A=10 -6 
K the process is largely non-adiabatic and the LZ mech- 
anism imposes coherent oscillations during the spin flip. 
They end in a Larmor spin precession around the z-axis 
as imposed by the longitudinal spin bias s, followed by 
a phonon-bottleneck and avalanche, as described above 
(fig. 2a). 

The same mechanisms are demonstrated by the con- 
tour plot in Fig. 3. Here, A is varied continuously be- 
tween 2-10 -7 K and 2-10 -6 K. Similar half-cycles as in 
Fig. 2a are obtained and shown in shades of gray. (S z ) 
varies from -1 (black) to +1 (white) and one can observe 
the sudden spin flip in zero field, followed by a phonon- 
bottleneck and avalanche sequence. The bottom half of 
Fig. 3 gives a zoom taken at very low and positive mag- 
netic fields. The spin precession is shown as alternating 
shadings; they gradually vanish at higher fields due to the 
limited coherence T 2 time or at higher values of A where 
the spin flip enters into the adiabatic regime. If the field 
is swept at higher rates and at lower temperatures, the 



p-4 



Coherent spin rotation in the presence of a phonon-bottleneck effect 




o 2x10* 4xio* ex to' feiuy" ixio* 
Magnetic Field H (T) 



Fig. 3: Contour plots of half-cycles, with the field ramped from 
negative to positive values. The (S z ) increases gradually from 
-1 (black) to +1 (white). The zoom shows coherent oscillations 
induced by the LZ mechanism. 

phonon avalanche can be strongly reduced. Thus the lat- 
tice phonons and the spins are brought to temperatures 
much lower than the cryostat one (adiabatic cooling). Such 
behavior has been observed experimentally in the case of 
the V15 molecular magnet [9]. After several field cyclings 
and for a crystal thermally isolated from the setup, the 
sample temperature drops significantly and the spin re- 
laxation time is increased from seconds to hours. 

In experiments involving ultra-low temperatures and 
scanning magnetic fields, active cooling of a spin system is 
therefore possible for a desired amount of time. It is im- 
portant to note that active cooling experiments have been 
performed with superconducting quantum bits ( [20], [21]). 
Quenching the lattice at lower temperatures using a PB 
mechanism could provide an additional knob to improve 
the active cooling of qubits or of any other quasi-spin sys- 
tem. 

Conclusions. — We present a simple phenomenologi- 
cal method to simulate the effect of non-adiabatic Landau- 
Zener excitations on spin reversal driven by a phonon- 
bottleneck mechanism. To this purpose, Bloch equations 
are written in the eigenbasis of the effective spin Hamil- 
tonian, considered to be a two-level system at low tem- 
peratures. The relaxation term is given by the phonon- 



bottleneck mechanism. To the expense of calculus time, 
the method can be generalized to multi-level systems, 
where the notion of Bloch sphere does not apply but the 
density matrix formalism is still applicable. 

The authors acknowledge support from the Florida 
State University, National High Magnetic Field Labora- 
tory (grant IHRP-5059 and DMR-0654118), NSF (grant 
CAREER DMR-0645408), DARPA (grant HR0011-07-1- 
0031 ) and the Alfred P. Sloan Foundation. 

REFERENCES 

[1] Wernsdorfer W., Adv. Chem. Phys. 118 (2001) 99. 

[2] Nellutla S. et al, Phys. Rev. Lett., 99 (2007) 137601; 

Bertaina S. et al, Phys. Rev. Lett. 102 (2009) 050501. 
[3] Thomas L. et al, Nature, 383 (1996) 145; Friedman J. R., 

Sarachik M.P., Tejada J., Ziolo R., Phys. Rev. Lett, 76 

(1996) 3830; Sessoli R., Gatteschi D.,Caneschi A., Novak 

M.A., Nature, 365 (1993) 141; Chiorescu I. et al, Phys. 

Rev. Lett., 85 (2000) 4807; Wernsdorfer W., Sessoli R., 

Gatteschi D., EuroPhys. Lett, 47 (1999) 254. 
[4] Wernsdorfer W. and Sessoli R., Science, 284 (1999) 133; 

Wernsdorfer W. et al, Physica B, 284-288 (2000) 1231. 
[5] Wernsdorfer W. et al, EuroPhys. Lett, 50 (2000) 552. 
[6] Chiorescu I. et al, Phys. Rev. Lett, 84 (2000) 3454. 

Chiorescu I. et al, J. Appl. Phys., 87 (2000) 5496; 

Chiorescu I. et al. J. Magn. Magn. Mat., 221 (2000) 103. 
[7] Waldmann O. et al, Phys. Rev. Lett, 89 (2002) 246401. 
[8] Rousochatzakis I. et al, Phys. Rev. Lett., 94 (2005) 

147204. 

[9] Chiorescu I. et al, Phys. Rev. B, 67 (2003) 020402(R). 
[10] Chen L. et al, Appl. Phys. Lett, 89 (2006) 252502. 
[11] Leuenberger M. N. and Loss D., Nature, 410 (2001) 789. 
[12] Bertaina S. et al, Nature 453 (2008) 203. 
[13] Abragam A. and Bleaney B., Electronic Paramagnetic 

Resonance of Transion Ions, (Clarendon Press, Oxford) 

1970, Chap. 10, pp. 541-583. 
[14] Saito K. and Miyashita S., J. Phys. Soc. Jpn., 70 (2001) 

3385; Rousochatzakis I. and Luban M., Phy. Rev. B, 72 

(2005) 134424; Vogelsberger M. and Garanin D.A., Phy. 

Rev. B, 73 (2006) 092412. 
[15] de Loubens G. et al, Europhys. Lett, 83 (2008) 37006; 

Chudnovsky E. M., Garanin D.A., Schilling R., Phys. Rev. 

B, 72 (2005) 094426. 
[16] Zener C, Proc. R. Soc. London, A137 (1932) 696; Lan- 
dau L.D., Phys. Z. Sowjetunion, 2 (1932) 46. 
[17] Van Vleck J.H., Phys. Rev., 59 (1941) 730. Stevens 

K.W.H., Rep. Prog. Phys., 30 (1967) 189. 
[18] Hellwege K.H., Koetzler J., Weber G., Z. Physik, 217 

(1968) 373. 

[19] Faughnan B. W., Strandberg M.W.P., J. Phys. Chem. 

Solids, 19 (1961) 155. 
[20] Valenzuela S.O.ei al, Science, 314, (2006) 1589. 
[21] Chiorescu L, Science, 314 (2006) 1549. 



p-5 



